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^^ . We investigate the origin of the kink effect (KE) in the relativistic 

Q ■ 

mean field theory by transforming the single-particle Dirac equation into a 

Schrodinger-like equation. It is found that relativistic self-consistent effects 

y . as well as contributions coming from the p meson determine the actual 

C ■ structure of the KE. However, the spin-orbit force generated by the p meson 



has no significant influence on the KE. 
PACS numbers: 21.30.+y, 21.10.Ft, 21.60.Jz 



The charge radii (r^) of the Pb isotopes have been measured a few years ago with 
a very high precision |]l|, their anomalous kink behavior being the most essential feature 
of these measurements. The kink effect (KE) means that the experimental data on Tc 
as a function of A in the isotopic chain display a change of the slope at N=126 with a 
gradual addition of neutrons. This phenomenon has been a matter of a detailed discussion 
by different theoretical groups during the recent years (see Refs. 01-0 and references 
therein). 

At the first stage, it has been shown that the density-dependent Hartree-Fock model 
with standard parametrizations of the Skyrme (SHF) [Q or Gogny [^ forces fails to 
reproduce the empirical charge isotopic shifts in the chain of Pb isotopes. 

On the other hand, in Refs. P,Q the anomalous behavior of the charge radii of these 
isotopes has been studied in the relativistic mean field approach (RMFA). The Pb data 
are well described by this type of theory. In Ref. IQ it is supposed that the success of the 
RMFA in reproducing the KE is achieved due to the weak isovector dependence of the 
spin-orbit force generated by the RMFA. 

In Ref. [0 the problem of isotopic shifts has been investigated with exhaustive 
detail, both in the SHF and RMFA. The authors of Ref. have proposed a new Skyrme 
functional SkI4 with a more general structure of the spin-orbit force than that of the 
standard Skyrme functional. They have also succeeded to reproduce the KE, but at the 
price of introducing an extra parameter for the two body spin-orbit potential. 

As the authors of Refs. 0,^ state, the standard RMFA provides an excellent de- 
scription of the anomalous kink in the isotopic shifts about ^osp^ without using extra 
fitting parameters. In the present paper we shall focus our interest into the problem of 
the origin of the KE in the RMFA paying special attention to the role of the isovector p 
meson in the kink structure. 

In the RMFA 0, the relativistic single-nucleon wave function ipl-f), for a nucleon 
with the rest mass M, satisfies the Dirac equation 

[-tha ■ V + (3{M + S{f)) + V{r)]ij{r} = S^(f), (1) 

where E = M + e and e is the binding energy eigenvalue while S{r) and V{r) are the 
scalar and vector potentials, respectively, defined as follows: 

S{r) = g^aif), 
V{r) = K.(r) + ToV.ir) + ^^K(r) = g^M^^ + r.g^Poi^ + e^^Aoir^. (2) 



In Eq. (2) a is the scalar meson field, ujq and po are the time components of the 
oj and p meson fields, respectively, while Vc is the Coulomb potential, g^, g^^ and gp are 
the coupling constants associated with the respective mesons. All meson fields satisfy the 
corresponding Klein-Gordon equations (see 0, for example). 

Eq. (1) can be reduced to a Schrodinger-like equation |^ which does not contain 
the first derivative of the wave function [0(r)]. It reads: 

[-^V^ + V;en*(r) + Vso{f^]m = e{l + ^)^(^1, (3) 

where the central potential Vcent and the spin-orbit potential Vso are given by: 

5-2 _ 't/2 y 

Vcentif) = 5 + r + ^^ + ^M + ^^-"*' (^) 



2M + £ + 5-\/' 



h^ 2. 



V5o(r) = — -W^/ ■ s. (5) 

IM r 

By writing the Dirac equation in this form, we are able to study separately the 
influence of the different terms entering Eq. (4) on the KE (we can analyze, for example, 
if their respective effects are of relativistic origin or not). 

Our calculations have been carried out for the Pb chain of isotopes around the ^^^Pb 
nucleus in the framework of the standard RMFA, with a linear model (L) |^ and two 
nonlinear models NL-SH [|10| and NL3 [|Tl| containing self-interactions of the scalar field. 



For simplicity, we shall mainly concentrate our discussion in the three Pb isotopes ^°^P6, 
^°^P5 and '^^^Pb. The most stable configuration of the "^^^Pb isotope in the two nonlinear 
models considered contains two neutrons in the liii/2 states. In the L model, the two 
configurations in which two neutrons fill up the l2ii/2 or 2(^9/2 single-particle states are 
almost degenerate (we shall refer to these two configurations as IC and GC, respectively). 
We approach the ^^^Pb ground state by the IC in both linear and nonlinear models, 
although the BCS formalism would be more appropriate in the first case. Nevertheless, 
we expect that our conclusions are also valid for the L model, at least, at a qualitative 
level. 

The self-consistent results for the charge radii corresponding to the exact models are 
summarized in Table I, case a, and in Fig. 1 (the r^ are normalized in this figure to the 



experimental value of the 208p^ nucleus). The KE appears in the three models considered 
but it is more strongly pronounced in models with smaller K modulus values. We notice 
that the KE is quite small for the GC, in spite of the fact that the 2(^9/2 state is less bound 
than the liii/2 state and that the radius of the 2(^9/2 state is larger than the radius of the 
liu/2 one. Thus, it is not the radius but the different neutron density associated with 
these two states that is essential to explain the different KE corresponding to IC and GC. 
The KE in the chain of Pb isotopes can be only explained if single-particle configurations 
including liii/2 states significantly contribute to the total nuclear ground state for the 
isotopes with A > 208. 

To understand how the kink is generated by the model we have begun from studying 
the influence of the different mesons on the KE. We have calculated the r.m.s. radii of 
different components of the single-particle potential: S, VL and Vp as functions of the 
number A, and we have established that it is only the value of Vp, i.e., the r.m.s. radius of 
the p meson potential Vp, that shows a pronounced kink while the other components do 
not. Since the p meson has an isovector character, we can expect that it plays a relevant 
role in the behavior of protons when neutrons are added to the nucleus. Thus, the result 
mentioned above confirms the relevant role of the p meson in generating the KE. 

Let us underline that the p meson contribution reveals itself in Eq. (4) in different 
terms of the total potential Vcenti^) + Vso{£) through the V component. Since the couphng 
constant Qp of the p meson is comparatively small, the most essential contribution of the 
p meson to the total potential comes from the (S+V) and (S*^ — V^)/2M terms, where 
the a and uj contributions cancel each other to a certain extent. 

To study the influence of the different components containing the p meson field on 
the charge radii of the nuclei in more detail, it is useful to make the following replacements 
in eqs. (4) and (5) corresponding to Vcentif) and Vsoif), respectively : Qp -^ XiQp in the 
(S + V) term, Qp -^ X2gp in the {S^ — V^)/2M term, Qp -^ x^Qp in the j^V term, gp — > X4gp 
in the AVcent term, and Qp -^ x^Qp in the the Vso potential. After that, we allow the 
parameters Xi_5 to vary in a continuous and independent way in the interval [0-1]. We 
have considered 4 different cases corresponding to the 4 following combinations of the Xi 
parameters: a) Xi^^ = 1 (already discussed) b) Xi <^ 1, X2-5 = 1, c) Xi^3_5 = 1, X2 = 0, 
and d) Xi_5 = (equivalent to take Qp = 0). 

The results for these cases are collected in Table I, and can be identified by the 
corresponding letter a, b, c or d in the last column. We have considered the models L, 
NL-SH and NL3, and the IC for the ground state. Fig. 2 shows the r^ results for the 
NL-SH model normalized to the r^ of the '^^^Pb nucleus. The case a contains the total 



p contribution and manifests the KE in good agreement with the effect experimentally 
observed. In the case 6, in which the p contribution to S* + V is almost eliminated, the 
kink has completely changed its structure in relation to the case a, now the kink seems to 
be, somehow, the mirror image of the kink in case a. This fact makes evident the crucial 
role the p meson plays in the kink structure through the contribution of the {S + V) term. 

Table I shows that r^ slightly increases when the parameter xi decreases. This fact 
is related to the decreasing of the proton binding energy with the decreasing of the Xi 
value (the opposite is true for the neutron energy). However the qualitative change of 
the kink structure cannot be attributed to this modification of r^ (which has only a small 
effect on the kink). Notice also that xi is taken equal to 0.2 and 0.1 (rather than 0) in 
the L and NL-SH models, respectively, to maintain the 3si/2 proton state bound enough. 

The effect of the p meson on the kink in the region 208 < A < 210 through the 
term S + V \n the potential is quite small for the more stable IC However, for the GC, 
a decreasing of the Xi parameter produces a sizable turn at the left side of this region of 
the kink. 

With regard to the case c in Fig. 2, the effect displayed is opposite to that shown in 
the case h and much smaller in magnitude. This result is a consequence of the opposite 
role the p meson field plays in the terms {S + V) and (S*^ — V'^)/2M of the central 
potential. For the GC the effect of the p meson on the kink through these two terms 
remains opposite to the case 6, but comparable in magnitude. 

Finally, for the case d, i.e., the case without p meson contribution, the behavior of Tc 
is intermediate between that corresponding to the cases a and h as it can be appreciated 
in Table I and Fig. 2. 

We have verified that the p meson contribution through the Accent and Vso terms 
has only minor infiuence on the kink. As we have explained above, this behavior is 
expected because the relative contributions of the p meson to the {S — V) term, entering 
Vcent and Vsoi are very small. 

To understand better the nature of the KE it would be useful to examine the partial 
contributions Tc^ of different single-particle states " /c" to the total charge radii r^. of the 



nuclei under consideration (re = Jrl^ + ... + r^^). Fig. 3 displays these values for the six 
bound states (closest to the Fermi level) calculated within the NL-SH model for different 
nuclei in the region 200 < A < 214. We see that the total value of r^ for A > 208 is 
mostly determined by the contribution of the l/iii/2, 15*9/2 and 1^7/2 proton states, rather 
than by that of the 2^3/2, 2^5/2 or 3si/2 proton states. Let us mention also that in the 
nuclei considered there are more protons of the first type than those of the second type. 



Moreover, the values r^k, for the states of the first type, are increasing functions of A 
for A > 208 while the respective values for the states of the second type almost remain 
constant with the increasing of A in this region. One can effectively describe this type 
of behavior as a result of the repulsive interaction experienced by the l/iii/2, ^99/2 and 
^97/2 protons which pushes them outside the central part of the nucleus and increases the 
value of Tc- From Fig. 3 it is also seen that the behavior of the total value of r^ for the 
nuclei in the interval 206 < A < 208 is strongly affected by the Vck value for k = 3si/2. In 
fact, a similar influence is also shown by the 2si/2 state, in spite of its very large binding 
energy. In cases where the p contribution to the {S + V) term in the central potential is 
reduced (cases b and d in Fig. 2), the contributions of the two 2si/2 and 3si/2 states to 
Ar~ = Tci^^^Ph) — Tci^^'^Ph) are even larger. The two neutrons added to the 3pi/2 state 
appear to be very effective to pull on the protons in these two states. 

Table I and Fig. 2 also show that in the region 206 < A < 208 it is easier to polarize 
protons by adding neutrons when protons are less bound (cases h and d), in accordance 
with the reasons given in Ref. |P|. However, the behavior of re in the region 208 < A < 210 
is the opposite. The r^ behavior shows that the cases that produce larger variation of 
Ar^ are less effective contributing to Ar+ = VcC^^^Pb) — Tci^^^Ph). This is a reflection of 
the fact that it is more difficult to increase the r^ value by adding extra neutrons when 
the radius is larger than in the case when it is smaller and when the neutrons previously 
added are more effective to increase Tc- 

As explained above, although the p meson contribution to the total energy is rather 
small in comparison with that of the a and uo mesons, it is just the p meson that is 
mainly responsible in the RMFA for the kink structure. However, even without p meson 
contribution (gfp = 0) there is a residual kink, case d in Fig. 2. In this case, we have 
investigated the role of the different terms entering Vcent and of Vso in producing the 
residual kink. To do that, we have artificially changed the contribution of these terms in 
a similar way as we did it for the different components of the p meson, but readjusting 
now the strength of the {S + V) term, if necessary, to keep the Tc values close to the 
experimental ones for the Ph isotopes studied, (the same conditions are used for the three 
isotopes). These studies reveal that the terms e-^, AVcent and Vso contribute in a slight 
and positive way to the kink structure, whereas the term (5*^ — V'^)/2M has an opposite 
effect that approximately cancels the effects mentioned above. We notice that taking 
about the 70% of the strength of the (5* + V) term entering Eq. (3) as the total potential, 
the KE (obtained in the Pb isotopes studied taking 9p = 0) is mainly reproduced as 
an effect of the self-consistent procedure. For example, if the S and V potentials are 
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approached by the Woods-Saxon functions with the radius proportional to A^^^ the KE 
disappears, although different diffuseness parameters for the potentials are used. 

We have studied above the role that Vcent and Vso play in generating and forming 
the structure of the kink. Another crucial relativistic effect comes from the fact that the 
source term in the scalar field equation, the scalar density (ps) , is smaller than the actual 
baryon density (pb). This difference of two densities is just the mechanism that brings 
about the saturation in the RMFA, and essentially distinguishes relativistic models from 
nonrelativistic ones. Thus it would be interesting to check also the role of this mechanism 
in the KE. One naive way would be to replace ps in the scalar field equation by ps- Since 
the nucleus would become unstable in this case, we add to Vcent in the single-particle 
Schrodinger equation (1) a repulsive density-dependent term formally identical to the 
term proportional to t^ appearing in the corresponding equation for the Skyrme-Hartree- 
Fock model. The coefficient t^ is chosen to fit the Tci^^^Ph). Under these conditions we 
observe that in two cases, with and without the p meson (cases a and d, respectively), 
the quantities r^ and r^ are reduced and so is the ratio t^ /t^ ■ Thus, the KE is reduced 
about 15% in the case a and 60% in the case d. The residual KE without the p meson 
almost disappears. 

Our conclusions can be summarized in the following points: 

1. The KE in Pb isotopes chain appears to be a general feature of the RMFA. It 
is observed in both linear and nonlinear models, being more strongly manifested in the 
models with smaller values of K. We have found that the p meson plays an essential role 
in forming the KE. In the case without the p meson a small kink remains that can be 
attributed to relativistic as well as to self-consistent effects. 

2. The KE mainly comes from the (partially) destructive interference of contribu- 
tions of the p meson to the {S + V) and (S*^ — V'^)/2M components of the single-particle 
potential. The {S + V) term leads to a KE of positive sign while the {S"^ — V'^)/2M compo- 
nent produces a KE of negative sign. The total effect has a positive sign and corresponds 
to the trend observed empirically (see cases h and c in Fig. 2). 

3. Excluding the p contribution from the two components of Vcent mentioned in 
the previous point, i.e., taking Xi = a;2 = 0, one gets a strong reduction of the KE. This 
situation is quite similar to that in which the p meson contribution is completely excluded 
((^^=0, case d in Fig. 2). The residual KE in this case is related to the self-consistent 
procedure. This effect completely disappears if one uses a Woods-Saxon potential in the 
Dirac equation rather than a self-consistent calculation. 

4. The contribution of the p meson to Vso has no influence on the KE. 



5. The relativistic effect that produces saturation in the RMFA provided by the 
difference between the scalar density and the baryon density, also contributes to a certain 
increase of the KE with respect to the nonrelativistic case in which saturation is produced 
by a repulsive density-dependent term as in the Skyrme-Hartree-Fock model. 

In a forthcoming publication we shall investigate the problem of the KE in the 
relativistic Hartree-Fock framework. 

One of the authors (L.N.S.) is very grateful to the University of Cantabria and to 
the Technical University of Munich for hospitality. This work has been supported in part 
by the DGESIC grant PB97-0360. 
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FIGURE CAPTIONS 



Fig. 1. The charge radii of Pb isotopes obtained with the L 0, NL-SH |T^ and NL3 
TT| models, normahzed to the experimental value of the ^osp^ nucleus. The empirical 



values (EXP) are taking from Ref. [^ 



Fig. 2. The same as in Fig. 1 for the NL-SH |TIj] model, and the four cases considered in 
Table I: a: a;i_5 = 1 (the exact model); b: Xi = 0.1, 0:2-5 = 1; c: Xi^3_5 = 1, X2 = 0; d: 
Xi_5 = (without the p meson contribution). 

Fig. 3. The single-particle r.m.s. radii for the proton states, as indicated, in the Pb 
isotopes obtained with the NL-SH |T0[ model. 
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TABLES 



Model 


rc(fm) 


Ar+(fm) 


Ar-(fm) 


Ar+/Ar^ 


case 


L 


5.455 


0.016 


0.010 


1.59 


a 


NL-SH 


5.502 


0.020 


0.010 


1.93 


a 


NL3 


5.517 


0.023 


0.010 


2.37 


a 


L 


5.577 


0.016 


0.018 


0.88 


b 


NL-SH 


5.651 


0.018 


0.022 


0.82 


b 


NL3 


5.718 


0.019 


0.022 


0.85 


b 


L 


5.421 


0.016 


0.008 


1.99 


c 


NL-SH 


5.477 


0.020 


0.008 


2.46 


c 


NL3 


5.495 


0.024 


0.008 


3.08 


c 


L 


5.510 


0.016 


0.014 


1.17 


d 


NL-SH 


5.558 


0.018 


0.014 


1.28 


d 


NL3 


5.568 


0.021 


0.013 


1.58 


d 


Exp 


5.503 


0.019 


0.011 


1.73 





TABLE L The charge radius (r^), Ar+ = rc{^^^Pb)-r^{^^^Ph), Ar" = rc{^^^Pb)-rc{^^^Pb) 
and the ratio Ar+/Ar^ for the hnear model L [9] and the nonlinear models NL-SH [10] and NL3 
[11], for the cases: a: zi_5 = 1; b: xi = 0.2 (L), 0.1 (NL-SH), (NL3), X2-5 = 1; c: xi,3_5 = 1, 
X2 = 0; d: xi_5 = 0. The last line contains the experimental results [1]. 
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